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The present status of the (JR) dynamical parton distribution functions is reported. Dif- 
ferent theoretical improvements, including the determination of the strange sea input dis- 
tribution, the treatment of correlated errors and the inclusion of alternative data sets, are 
discussed. Highlights in the ongoing developments as well as (very) preliminary results in 
the determination of the strong coupling constant are presented. 

The dynamical parton distributions of the nucleon at Q 2 > 1 GeV 2 are QCD radiatively 
generated from valencelik^\ positive definite input distributions at an optimally determined 
low input scale Q 2 , < 1 GeV 2 . Thus the steep small- Bjorken-x behavior of structure functions, 
and consequently of the gluon and sea distributions, appears within the dynamical (radiative) 
approach mainly as a consequence of QCD-dynamics at x < 10~ 2 [1 . Alternatively, in the 
common "standard" approach the input scale is fixed at some arbitrarily chosen Qq > 1 GeV 2 , 
and the corresponding input distributions are less restricted; for example, the mentioned steep 
small- a; behavior has to be fitted. 

Following the radiative approach, the well-known LO/NLO GRV98 dynamical parton dis- 
tribution functions of [2] have been updated in [3], and the analysis extended to the NNLO 
of perturbative QCD in [4]. In addition, in |3j 0] a series of "standard" fits were produced in 
(for the rest) exactly the same conditions as their dynamical counterparts. This allows us to 
compare the features of both approaches and to test the the dependence in model assumptions. 
The associated uncertainties encountered in the determination of the parton distributions turn 
out, as expected, to be larger in the "standard" case, particularly in the small- a; region, than 
in the more restricted dynamical radiative approach where, moreover, the "evolution distance" 
(starting at Q\ < 1 GeV 2 ) is sizably larger Oil]. 

The strong coupling constant a s (M§) was determined in our analyses together with the 
parton distributions, in particular it is closely related to the gluon distribution which drives 
the QCD evolution and consequently its uncertainty is also smaller in the dynamical case. We 
obtained a s {M%) = 0.1124 ± 0.0020 at NNLO, and 0.1145 ± 0.0018 at NLO in the dynamical 
case; to be compared with a s (Mf) = 0.1158 ± 0.0035 at NNLO, and 0.1178 ± 0.0021 at 
NLO in the "standard" one. The difference between these values is to be interpreted as a 
genuine uncertainty stemming from parameterization dependence, e.g. our contribution to the 
next PDG determination will be an average over NNLO dynamical and standard results. We 
consider this approach to be more realistic than considering only the reduced errors stemming 
exclusively from experimental uncertainties in each single analysis [5] . 

1 Valencelike refers to aj >0 for all input distributions xf(x,Q^) ac x a f(l — x] f, i.e., not only the valence 
but also the sea and gluon input densities vanish at small x. 
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The dynamical predictions for Fl(x,Q 2 ) become perturbatively stable already at Q 2 = 2 - 
3 GeV 2 , where precision measurements could even delineate NNLO effects in the very small-x 
region. Moreover they are positive and in excellent agreement with the latest HI data [BJ. This 
is in contrast to the results based in the common "standard" approach, which are less precise 
and in some cases turn negative at the lower Q 2 values. 

The inclusive production cross sections for W + , W~ and Z°-bosons form important bench- 
marks for the physics at hadron colliders. Our NLO [7] and NNLO [5] predictions for these 
standard candles processes were presented not long ago, and a detailed comparison of the predic- 
tions predictions based on the available NNLO parton parameterizations was performed in [9]. 
According to these studies the rates for gauge boson production at the LHC can be predicted 
with an accuracy of better than about 10% at NNLO. NNLO predictions for the Higgs boson 
production cross sections for Tevatron and LHC energies were also considered; the production 
rates could be predicted with an accuracy of about 10-20% at the LHC. The inclusion of the 
NNLO contributions is mandatory for achieving such accuracies since the total uncertainties 
are substantially larger at NLO. 

During the time passed since the last determination of our PDFs there have been experi- 
mental and theoretical developments which deserve consideration. From the experimental side, 
besides results coming from LHC which we will not discuss here, the combined HERA data 
on neutral current and charge current DIS has been published [10] ; the combination of the 
charm production data is ready but unfortunately still not publicly available. For the older 
fixed-target data in our analyses we use now directly the measured cross-sections (instead of 
the extracted structure functions) since potential problems with some of these extractions have 
been pointed out |llj . A wealth of deuteron data have also been included, for which description 
we use the nuclear corrections of |12j . Further theoretical improvements include the use of the 
running mass definition for DIS charm and bottom production 13 . 

In this proceeding we will concentrate in several concrete aspects of the ongoing update, 
modifications in the ansatz for the strangeness input distributions, and some consequences of 
the treatment of correlated systematic errors, in particular of the normalization uncertainties. 
Some preliminary results on determinations of a s and the (future) inclusion of higher- twist 
terms in our analyses are also briefly discussed. 

As in previous dynamical determinations |2] , since the data sets used are insensitive to 
the specific choice of the strange quark distributions, the strange densities of the dynamical 
distributions in [3j H] have been generated entirely radiatively starting from vanishing strange 
input distributions: 

s(x,Q 2 ) = s(x,Q 2 )=0 (1) 

at the low input scalcQ The plausibility of these assumptions was investigated in |14j by 
confronting the predictions derived from dynamical distributions determined in this way with 
the most precise dimuon production data |15j , which were the most sensitive to the strangeness 
content of the nucleon until the recent ATLAS measurements [16] . A good agreement was found 
at NLO. Since the NNLO corrections (thus the experimental acceptance corrections) for these 
process are not known, a consistent study of these data at NNLO is not possible. However 
it is apparent that unless the NNLO corrections were very large (which seems unlikely) our 
JR09 NNLO distributions would generally undershoot the data. Indeed the would-be NNLO 
results favor s ~ 0.1 (u + oTj at the input scale. In order to estimate the implications which this 

2 In the "standard" case, where Q 2 > 1 GeV 2 , the strange input distributions were chosen s(x, Qfj) = 
s(x,Qg) = | (u(x,Qq) + d(x, Qq)) > as ' s conventional [3][3]. 
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could have for our predictions for benchmark processes at the LHC we have repeated our JR09 
dynamical analysis using s(x,Qq) = s(x,Qq) — \ (u(x, Qq) + d(x, Q 2 ,)) as in our "standard" 
fits. The resulting input distributions are compared with the original JR09 input in Fig. [I] It 
can be observed that the increase in the strange distributions is compensated by a decrease in the 
light sea, while the other distributions remain practically unaltered. The mentioned benchmark 
cross-sections for W and Z production are practically identical in both cases. Nevertheless the 
determinations of strange input distributions will have to be revisited in light of the mentioned 
ATLAS measurements 116] . which seem to imply a greater strange content of the nucleon. 

Another aspect of our analyses in which we have been working is a more complete treatment 
of the correlated systematic uncertainties of the data. The least-squares estimator that we use 
to take them into account has been explicitly written down in Appendix B of |17| . This is used 
as well for normalization uncertainties, and implies that data and theory are relatively shifted 
by amounts which are determined by the minimization but not restricted further. This is in 
contrast to the treatment that we used in our (G)JR analyses, where the normalization shifts 
were limited by the experimental normalization uncertainties, and has resulted in small shifts 
on the valence distributions in the medium x region. 

Yet another issue raised by the ABM collaboration is the necessity of including higher twist 
for the description of fixed target data (18j : even if moderate kinematic cuts are used to select 
the data included in the fits; in particular for SLAC and NMC data. The kinematic cuts in 
our G(JR) analyses were Q 2 > 4 GeV 2 and W 2 > 10 GeV 2 , and were applied to the F 2 values 
extracted from different beam energies and combined. The description was good for NMC data 
and rather poor for SLAC data. However since the number of points for these experiments was 
rather small (about 100 data points for NMC and 50 for SLAC), the values of the cuts did not 
affect very much the results of the fits. 
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This picture changes if, as the ABM collaboration does, data on the cross sections for 
individual energies are used, which amounts to hundreds of data points for each experiment. 
As mentioned in the introduction, we also use these data in our current preliminary analyses, 
for which the virtuality cut has been raised to Q 2 > 9 GeV 2 . We find results which are rather 
similar to our JR09 analyses, with a s (M§) values about 0.113 to 0.114, depending on the input 
scale. If the usual cut of Q 2 > 4 GeV 2 is applied to these data a s (M§) raises to 0.1176; a 
tendency in agreement with the ABM observations [18]. The inclusion of higher twist terms 
in the theoretical description should strongly reduce the dependence of the outcome of the fits 
and is currently under consideration. 

To summarize, we have discussed several aspects of the ongoing update of the JR dynamical 
parton distributions, in particular the determination of strangeness input distributions, the 
treatment of normalizations, and results on a s (Af§) and their dependence on kinematic cuts 
and the treatment of higher-twist contributions. Changes in the parton distributions induced 
by all this modifications in the analysis and in the data have been, until now, quite modest. 
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